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ABSTRACT 

The objective of the present investigation was 
to arrive at an optingil design of the basic structuial 
element namely the Tee beam. The limit state theory, 
a recent design method was adopted. The optimization 
was carried out for the minimum cost . A nonlinear 
programming method with interior penalty function was 
utilized . A great deal of emphasis was given to the 
serviceability conditions like deflection and crack 
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CHAPTER I 
IITTRODUCTION 


1 . 1 GENERll 

The "basic role of a structural engineer is to plan, 
design and construct a structure. A structural designer 
in particular should "be equiped with the loiowledge of loads, 
materials, design methodologies and the advances in compu- 
tational techniques. Traditionally structural designs were 
evolved over a period of time. Specially during the II 
world war, the acute shortage of materials like steel made 
the designers in the war tom countries to take a fresh 
look at the reserve strength of the materials. This has 
led to the plastic design of structures. However, the 
advent of computers has revolutionized the hitherto unsolved 
pro"blems of the early fifties and la,rge structures can "be 
designed effectively for minimum cost. 

1.2 DESIC2T PROCESSES 

The design of a structure means the proportioning 
of various components according to the needs of the customer, 
A creative sense, scientific knowledge and a little 
practical esperience when combined, produce a satisfactory 
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solution. A thorough knOY/ledge of the properties of 
materials are as important as the methods of structural 
analysis. Different idealizations of the material 
properties describe different design methods. The object 
of structural design is to obtain a structural solution 
which gives the greatest overall economy by providing the 
optimum while satisfying all the requirements of the 
stmcture. It is clear that the structural costs rise 
rapidly if other requirements become too onerous. So, 
a method v;hich is simple, rational aiid readily applicable 
should be the aim of a designer. To understand the general 
design tools available at present, a brief description is 
given in the following paragraphs, 

1.2.1 ELiSTlC DESI&h (1) 

This is the most widely used traditional and well 
established design tool. It assimies the elastic behaviour 
of the stiucture associated with the assumption that the 
materials used individually and as a composite usually obey 
the Hooke's lav/. A structure is generally designed for a 
particular case as flexure and other aspects are checked 
so that the design is feasible. In this method, the 
factor of safety is ensured by limiting the maximum stresses 
under working loads to a fraction of the values determined 
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by experiments. Composite materials are transformed as 
unitary materials by using the modular ratio of the materials. 
The section may thus be transformed into a homogeneous one for 
the purposes of analysis. Unfortunately, the actual factor 
of safety i,e, the ratio between the collapse load and the 
working load, is different for different structures and 
loading conditions. Hence this design method is not 
rational so far as the actual factor of safety is concerned. 
However, this method is handy for studying the behaviour 
of structure under working loads. 

1.2.2 ULTIMATE LOAD DESIGF 

This method is rational when compared with the 
elastic method described earlier. In this method, the 
ultimate values such as ultimate moment of capacity-^, 
ultimate shear capacity etc, are computed. The stress- 
strain curves foi? concrete and steel are idealized prior 
to their application in the design process. For example, 
the actual stress block in compression zone of a concrete 
section is idealized as a rectangular stress block although 
the actual stress-strain diagram is rectangular parabola 
which inturn is a better approximation of the recorded 
stress-strain curve of concrete in compression. The 
ultimate strain in concrete is assumed to be 0,003 for 
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all grades. The steel stress-strain curve is idealized as 
a bilinear curve. The slope suddenly changes to zero at 
the yield point. In this method, load factors axe used to 
arrive at an ultiraate load from the given working loads. 

Various types of loads (dead load, live load, wind load 
etc.) have diffe.vont load factors. In the design of beams, 
usually the section is proportioned for flexure. This 
method, however does not give any idea about the behaviour 
of the structure at working loads. 

1.2.3 LIMIT STATE DESIGN 

The recent tool of design available to an engineer 
is the ’Limit state design’. It is more rational and an 
advanced design method wherein more accurate stress- strain 
curves are used in the design of the structure. The term 
limit state -is , ; defined in the f»i.b©utrig' Clucptcr, 

The advantage of this method is that the limit state can 
be choosen according to the needs of the customer. It is 
also fairly easy to apply this method for ordinary design 
problems, Eor the first time in literature a detailed 
design for the serviceability requirements like deflections, 
cracks, local damages etc, are considered in this method 
of design. Recently, codes have also started adopting 
this as the standard design. This method of design adopts both 
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the elastic and ultimate theories that woedescrihed earlier. 
More accurate determination of design loads is possible with 
the introduction of the partial safety factors, iULl this 
has made the author to adopt this process of design in the 
present thesis, ilthough this is byfar the best design 
method, it does not yield unique solutions. 

1.2.4 PROBiBIIISTIC DESIGN 

In nature, the material properties and the environ- 
mental forces are not deterministic. But in the foregoing 
design methods, the loads and material properties are 
considered as doteministic. The purpose of any design is 
the achievement of the probabilities that the structure is 
fit for use under all the circumstances and loading condi- . 
tions. It is evident that any general design requires some 
quantities of the design variables obtained by experimen- 
tation, They are found out applying the statistical methods. 
But in practice the strength varies. Such a design incor- 
porating the variability of materials and loads • , ' ;■ 

' may be termed as ’Probabilistic design’. This 
branch of science has caught the attention of many researchers 
around the globe. The basic design process hinges on the 
decision of the engineer as to how much level of probability 
is to be allowed for in the design. The general design 
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equation looks as follows. 

P [f (5) I < p (1.1) 

which is read as the prohahility (p) of a design function 
(f(D)) is at most limited to the probability level (p) , 

1.2.5 OPTIII£L DESIGN 

The designs that are arrived at are nonunique. Hence 
some design are bound to be uneconomical in terms of material 
consumption, cost etc. Thus it is imperative for optimizing 
a given design process in a systematic way. The optimal 
design may be in any domain of the material behaviour as 
elastic, plastic, visco-elastic etc. This is achieved either 
by using mathematical programming or optimal synthasis 
tools. 


1 . 3 OPTIMIZATION NEEDS 

In the case of aerospace vehicles minimum weight 
is the key factor whereas in developing countries like 
India, there shoifLd be a consceintious effort on the part 
of a designer to see that the structure consumes not only 
the minimum quantities of materials but also the overall the 
cost also should be minimum. The disproportionate increase 
in costs of building materials will make an engineer think 
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about the design adopted. Large structures (multi- storied 
buildings, dams etc.) can now be designed in an optimal 
manner, because of the advent of the fast digital computers, 

in optimal design problem has an objective function 
and certain constraints which are expressed as functions of 
design variables. 

1.4 LITERifURE SURVEY 

The elo.stic theory of stmctural design being 
conservative, the need for rational and less conservative 
methods of design was felt by structural engineers as 
early as 1937. The concepts of ultimate load theory were 
developed by assuming different material properties to suit 
their convenience. One of the researchers that had developed 
the theories, C.S. i^itney was the most successful one. 

Eor the first time he had proposed the design by the consi- 
deration of a stress block in the compressive zone of 
concrete* He had assumed a simplified shape of the stress 
block which was rectangular (1)* Later researchers had 
conducted experiments and different propositions for stress 
blocks were given. Luring that period engineers had not 
taken care the serviceability conditions. They had proposed 
designs for collapse of the stmcture. Researchers had 
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started thinking about the analysis for serviCvSability 
and local damages around 1965. 

Bengt 3, Broms (2) conducted experiments on short 
tension members and developed a simple method for calculation 
of crack width and crack spacing in R. 0. members. It was 
felt afterwards that a rigorous analysis method had to be 
propesed, unfortunately no researcher could give exact 
determination of crack widths as the nature of the cracks 
is irregular. Then people had started to control the 
cracks by same means. The imerican Concrete Institute 
Committee submitted a report (3) in 1972 on the control of 
cracks. It covered the basic mechanism of microcracking 
and fraction mechanics in concrete. This report had made 
a mention of the control of cracking due to shrinkage and 
flexure . It had also reported the long-tem effects on 
cracking. In this report a recommendation for construction 
was also provided in order to control the cracks. Gergely (4) 
had studied the effects of distribution of reinforcement 
for the control of crack. Later researchers Fawy (5), 
Orenstein (6) had suggested methods for control of cracks 
in two-way slabs. ELveny (7) had studied the behaviour of 
R.C. monber under normal conditions of drying to predict 
shrinkage stresses and deformations. Theoretical deter- 
mination of curvatures were suggested. Deflection of 
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R. G. flexural members was also studied by Branson (8) and 
principal factors that effect the short-teim and the long- 
term deflections were discussed. jO-Cock and Pauw (9) had 
proposed design methods for controlled- deflection design 
for R. C. beams and slabs. Keeping all this in view 
Hughes( 1 ') had given detailed description of the limit 
state theory and proposed certain design methods. MacGinley 
in his book (to) presented design theory and worked •. a 
number of examples based upon the British code (11) which 
had adopted the limit state design method for the first 
time. The Indian Standard Institution is also bringing out 
a code cf practice based upon the limit state theory (13), 

1.5 SCOPE OE THE THESIS 

In this investigation, the optimal design of Tee 
beam has been taJcen to develop ready reckoners for obtain- 
ing near optimal preliminary design for office floors. 
Definitions, of all of the terms used, are presented in 
Chapter 2, A brief mention of the method is also given in 
the same chapter. OHie design applied to a particular case 
of Tee beam is presented in Chapter 3. 0?he advanced method 
of design named the limit state method is applied. Hie 
optimization is done for cost to get economical designs of 
Tee beams. The Tee beam is selected as the element for 



1-10 


design v/ith a view that it is the most widely adopted 
structural el em jnt * 


The results and discussions are presented in 

(m. 

Chapter 5 andj^overview is presented in Chapter 6, 
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CHAPTER II 

DSEIwITIOHS AM) METHOPOIOGY 

2.1 GEHERAI 

Every scientific process has certain basic terms 
v/hich are to be defined, Eor the sake of con^jleteness 
definitions of the terms used in. the present work are 
given herein. It is also necessary to give a brief account 
of the methods applied, 

2.2 DSPIRITIONS 

2.2.1 STRUCTURE 

The skeleton of frame work, with the primary 
elements as beams, slabs, columns and walls, to carry the 
permanent and superimposed loads safely. It transmits 
the total load to the foundation, 

2.2.2 STRUCTURAl ELE’TMTS 

(1) Beams: Horizontal members carrying roof and floarr 
slabs. They resist loads in bending, shear and bond 
and may be simply supported, continuous or cantilevered. 

(2) Slabs: Horizontal plate elements vhich carry floor 
and roof loads. They may be simply supported or 


continous 
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2.2.5 S jPPOHT C OiTI) I2I OIS 

(1) Simply supported: Mechanical hinges which do not 
resist any raoment, transfering the loads to the 
supporting members, 

(2) Oontinous: Members which are continuous over the 
supports have these support conditions. 

( 3) Fixed: Ihe support conditions wherein the member 
is built into the support, in example of this is 
a Cantilever beam, 

2,2.4 MAlEBIiHS ilTO PROPERTIES 

(1) Bond: The adhesion of concrete to reinforcement, 

(2) Reinforced concrete: The concrete along with a rein- 
forcing material acts as a composite material. Used 
in the constructions of framed and walled buildings, 
bridges, culverts etc, 

(5) Creep: Strain or plastic flow that occurs in concrete 
due to a sustained load which is expressed as a 
creep coefficient (0) in case of concrete. 

(4) Cube strength: Crushing strength of a cube of 
150 mm sides, a.t 28 days. 

(5) Shrinkage: Contraction of concrete due to temperature 
changes. This is irreversible and results in cracking 
of the concrete. 
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(6) Strain(s): The change xn length per unit length of 
a member. It can be caused by physical influences 
shrinkage etc. 

(7) Stress- strain curve: A plot of stress against strain 
which shows the behaviour of the material under 
physical inflxiences (loads). 

(8) Yield Strength(fy.) : The limit of elastic region, 
beyond which there is a rapid increase in strain 
for a small incremental load and or for no increase 
in load. 

(9) Characteristic cube strength • The strengths of 

materials below which not more than 5 percent of test 
results fall. 

(10) Young’s modulus (E) : The ratio of stress to strain for 
a linear - elastic material. 

2.2.5 COMPOmTS or A TEE ESM 

(1) Tee beam: A rectangular beam which is cast monolithic 
with the slab over it. It shapes out into a ’T', 

(2) Elange: The portion of the slab whose width (B) is 
in effect to act monolithic ally with the rectangular 
rib below it. 

Rib: The rectangular beam lying underneath the slab, 

whose width is denoted as'b' . 


( 5 ) 
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(4) Depth (D) : The total depth of the rih and the 
slab put together is called the depth of the 
Tee beam, 

(5) Effective depth (d) : The depth measured from 
the top of flange to centre of gravity of the 
tension steel in the beam. 

(6) Effective cover: The amount of cover of concrete 
measured as the depth of concrete from the nearest 
face to the centre of gravity of the tension 
steel (Ag-p which is denoted as dg and as d^^ in 
Case of conpressive steel in the beam. 

(7) Span: The centre to centre distance between the 

supports of the slab (1„) , resting on the beams ■' 

or of the beam (l-j^) , resting on the walls or columns, 

2.2.6 LIMIT STATE METHOD 

(1) Limit state: State at which a structure becomes 
unfit for use. The main limit states are the 
limit states of collapse and serviceability. 

(2) Characteristic loads: The vyorking or service loads, 
classified into dead and imposed loads. They have 

a low probability of being exceeded during bhe 
life time of the structure. 

(3) Design loads: The loads obtained, from the characteri- 
stic loads, by multiplying the same by the partial 



2.5 


safety factors of safety for loads, 

(4) Design strength: Ihe characteristic strengths 
divided hy the partial safety factors of safety 
for materials, 

(5) Partial factors of safety (a): Factors applied to 
loads and materials to obtain the design values. 

They allow for uncertain! ties in the estimation of 
loads and variations in the strengths of materials. 
These control the safety of structures, 

(6) Ultimate strength of sections: A value based on the 
idealised stress-strain curves for the materials, 
taking account of the plastic behaviour. 

(7) Modular ratio (m) : The ratio of the Young’s modulus 
of steel to that of concrete, 

(8) Analysis: The process of determining the design 
actions, that are the shear, moment, bond etc, in 
the present work, 

(9) Moment of resistance( li^) : The moment that can be 
safely resisted by a structural element, 

(10) Deflection: The displacement of a member with 
respect to its original position. Generally measured 
with respect to the axis of the member. 

(11) Crack ’wLdth (¥) : The width of a crack fomed in 
concre1;e due to shrinkage, creep or fracture of 
the structure. 
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(12) Collapse: Ihe failure of a structural member, 

(13) Short-term, Reflections (y^^) : The deflection of a 
structure in a short duration. 

(14) Long-tem deflection ( 3 r_, v ) : The deflection of 
a structure v/hich are dependent on the time of 
loading. These are generally shrinkage and creep 
deflections, 

(15) Spalling of concrete: The ripping out of concrete 
due to excessive deflections there by causing 
cracks in concrete, 

(16) Dead loads: The loads that are permanent in nature 
contributed by the weight of the structural elanents, 

(17) Imposed loads; The loads lAhich are not permanent 
and sometimes may be dynamic in nature, 

(18) Moment of Inertia (l) : The second moment of area 
with respect to any axis. It is referred to an 
axis. 

(19) Stress block: The idealised stress distribution, 
of concrete in compfession, above the neutral axis, 

(20) Neutral axis: The axis where the strain in the 
cross-section is zero. 

(21) Depth of neutral axis (d^): The depth measured from 
the farthest fibre in compression to the neutral 


axis 
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(22) Under-reinforced sec cion: L section where the 
steel gets over stressed and therefore yields. 

The strain increases enormously without any increase 
in the stress therehy shifting the neutral axis 
towards the compression face. As a result the 
final failure of concrete occurs, which is a 
secondary effect. This failure is called a 
'Tensile failure' or 'Secondary compression 
failure ' , 

(23) Balanced section: The one for which the steel 
yields as the concrete just reaches its ultimate 
strain. That is, the beam fails simultaneously 
in flexural compression and tension. 

( 24) Over-reinforced section: The steel remains elastic 
while the concrete strains without increasing its 
moment. So, the neutral axis falls while the 
concrete i^apidly strains and crushing of concrete 
in compressive zone takes place. 

(25) Curvature (0) : Defined as the reciprocal of 
radius of curvature, ■which is the response of a 
structure due to the lateral loading that causes 
flexural stresses in the members, 

2.2.7 OPTIMIZATION 

(1) Optimisation: 3he process of improving a merit 
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function to get th.e oest solution to the given 
problem, 

(2) Design variables (x^^) : The quantities #iich are to 
be determined to satisfy the given norms of a 
design. 

(3) Constraints (g^) : The restrictions on the 
functions of the design variables which influence 
the physical nature of the problem, 

(4) Objective function (f) : The function that has to 
be extremized (minimized or maximized) to get the 
best and most economical solution to a given problm. 
Some times this is also called a cost function or 

a merit function. This can also be stated as the 
criterion with respect to which the design is 
optimized. This is expressed in terms of the 
design variables. 

2 . 3 IIITIIODOLOC-Y 

A summary of the different design methods available, 
for the design of a structure, are presented in the preceding 
chapter. The recnent method of design namely the 'Limit 
state design' is adopted in the present thesis, with 
facilities of the computer being available, the design is 
duly optimized for the total cost of the structure. 
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2.5.1 LIIHT STATE OP GOILASE 
2. 5.1.1 BEMDINC- 

The selected stress distributions axe detailed in 
tbe next chapter. The stress block factors K^,K2 and 
are determined as ( 1 ) . 

1:3 = (27 + 0.55 f^) / (22 4- f^) (2.1) 

Kg = 0,5 - 1^/550 (2.2) 

where is the cylinder ( 500 mm dia, 600 mm long) 
crushing strength. The relation between the characteristic 
cube strength and the cylinder crushing strength is given 
according the equation (5.3). 

Determination of the depth of neutral axis (d^) is 
esgplained in article 5.3.1.2(a). After computing the depth 
of neutral axis the rectangular parabolic stress distribution 
(Pig. 5,3c) is used to determine the compressive force. The 
total coir^jxessive force is divided into two parts and is used 
for coirputations of the moment capacity of a section, A 
detailed account of this is given in the following chapter. 

The bending moment (M) on the section by external 
loads is calculated as 

M = kg ^ 1^ (2.3) 

where k2 is a factor dependent on the support conditions 
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of the str^.ictiire. This thesis has incorporated a value 
of 0.1, is the ultimate load computed as 


u d d s s 


(2.4) 


where and are the partial safety factors to dead 
load and imposed load respectively. In the present work 
cc^ is taken as 1,4 and oc^ as 1,6 as specified by unified 
code of practice (11). 


2 . 3. 1 . 2 SHEAR 


The cliange in moment along the span of the beam 
Causes shear on the cross-section of the beam. This in 
association vd.th the bending stress induces diagonal cracks 
in the structure. Hence it is necessary to suitably design 
the structure taking care of the shear which in turn controls 
the diagonal cracking. The shear stress on the section is 
obtained in the same way as it is done in the case of an 
elastic method of design, except that the lever arm is 
replaced by the effective depth of the section, A section 
is also designed to see that the shear stresses caused or 
not too high and are limited to a maximum allowable 
shear stress (■'^ 3 j_]_) \’'^ich is related to the characteristic 

strength of the concrete cube as 


"^all “ ^cu 


(2.5) 
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Detailed analysis is given in the following chapter, 

2. 3. 1.3 301TD 

Bond is defined earlier as adhesion between the 
concrete and the reinforcing material. The bond stress 
analysis is done in the same way as it is done in the case 
of elastic theory. The bond stress varies along the length 
of reinforcement but for all practical considerations it 
is assumed to be constant. Bond in reinforcing concrete 
plays two different roles, Firstly, it provides the 
required anchorage to the reinforcement so that total 
bending stress is transfered to the reinforcement over a 
required length and secondly the average bond (v-^^) requires 
sufficient amount of perimeter of the reinforcing material 
such that there is no slip of the reinforcement. The 
latter stress is called the average bond stress. The Tee 
beam does not have the problem of anchorage length except 
where there is a negative moment and curtailment of the 
main reinforcement. 

The average bond stress allowed in concrete is 
dependent on the grade of concrete and the type of reinfor- 
cing bar used, TZhere plain bars are used without deforma- 
tions or ribs over it then the allowable average bond 
stress can be calculated by the expression, i ; 

Ci.'* •” '■ 

d'ik 
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^ba ^cu''^ ^cu N/mm^ 

= 0^04 + 0.30 for < 40 N/mm^ 

= 1.9 for > 40 N/mm^ (2.6) 

In case , deforned bars are used the®e values given by 
equation (2.6) should be increased by 60 percent. 

2.5.2 LIMIT- STATE OE SERYICSABILITY 

It is necessary that a structure should ser've 
the function for which it is constructed. The functional 
utility will be lost if the structure has irregularities 
like large deflections, wider cracks, vibrations under 
working loads etc. In such a case users experience psycho- 
logical fear for living in those dwellings. So, it is 
felt that a rigorous analysis of such secondary effects 
of the structure will enable an engineer to assess them 
properly and to limit them according to the needs of the 
customer. In analysis of the deflections and the cracks 
is done under this heading. The other serviceability 
conditions arc out of scope of this thesis. The partial 
safety factors each of 1,0 are used to get the design loads, 

2. 3. 2.1 DEPIEGTION 

The deflection of a reinforced concrete structure 
can be classified into two viz. short-term and long-term 
deflections. Tlie response of a structure immediately 
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after applice.tion of a load is short term deflection. 

In the other case any sustained load will make a concrete 
structure creep. It is also evident that concrete has 
shrinkage effects thereby strain the concrete further. So, 
these two effects of creep and shrinkage make a reinforced 
concrete structure to deflect and such a deflection is 
termed as long-term deflection. The procedure to analyse 
the structure for estimating the short and long-term 
deflections is presented in the following articles, 

2 . 3 . 2 . 1 ( a) SHORT- TERI'I DEPLEC TI ON 

As already mentioned reinforced concrete is not 
perfectly elastic. But for all practical conditions it can 
be safely assumed as an elastic material under the normal 
working loads. This assumption makes it simpler to analyse 
the structure easily and enables an engineer to apply the 
principles of elastic theory in the detemaination of 
deflections. 


The general relationship relating the deflection 
of a structure to the moment at any section is ( 1 ) . 


SI 


dx2 


- M 


( 2 . 7 ) 


vAiere y is the deflection of the strudture, 

M is the moment on the sectiem and 
p is the flextural rigidity of the structure 


considered, 
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But the Curvature of the sta:n;.ctural member can be 
approximated as 



On substitution and rearranging the terms one gets 

M 

0 ( 2 . 9 ) 

El 

The deflection y can be expressed, in terms of the 
curvature 0, and the length of the member 1 as 


7 =k 1^ 0 


( 2 . 10 ) 


v/here k = a constant given by the laws of mechanics 
applied to the structure. 

The constant (k) can b? easily calculated using the 
principles of mechanics or can be taken from standard 
tables which give its value for different boundary 
conditions. It is evident that one has to calculate the 
moment of inertia to compute the curvature and then the 
deflection of the structure. In the present analysis the 
moment of inertia ( Ir) of cracked section and the moment 
of inertia (ij of the gross section are calculated to 

o 

account for the stiffening effect of concrete in tension 
zone. The moment of inertia of the cracked section is 


Bd: 


I ” 


n 


-I- m 


( 2 . 11 ) 


3 
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whf’re the o.epth of neutral axis is computed from the 
following equation 


0.5 Bd^4- 


( 2 . 12 ) 

If the coimputed depth of neutral axis is greater than the 
depth of flange, then the depth of neutral axis is 
computed from equation (2.13) and correspondingly the 
moment of inertia of the section is obtained from equation 


(2.14). 


0.5 3a^+ 


|^( B-b) d^4- (m- 1) AgQ + 


( B-b) p 

= — d^ + (m-l)Agg d^ + (2.13) 

bd^ (B-b)d^ 'if 2 2 


+ 


12 






(2.U) 


The grcBs moment of inertia is calculated using equations 
(2.15) through (2.18) with the depth of neutral axis being 
computed from equation (2.15) in case the depth of neutral 
axis is greater than the depth of the flange; if not, it 
is coii 5 )utcd from the equation (2.17). The moment of inertia 
is calculated by one of the two equations (2.16) and (2.18). 
bd 

d = — 
n 2 


/ (Bd^-bd^ + bB) 


(2.15) 



2.16 




'§ 


12 


■— + (&-b)cL^ (d^ - 


>)2 


bd; 


+ 


n 


^ --.-(D-d„)' 


(2.16) 


djj = 0.5 (bD^ + Bcl| -bd|) /(bB + Bd^-bd^) 


^S = 


B^b ^(B-b) 

(]>-dj5 + -- — (d a )3 
3 f 


(2.17) 


3 n 


TV 


(2.18) 


Then the effective moment of inertia ( I_) which encounter^ 

t? 

the stiffening effect of concrete in tension zone is 
Calculated as 


^e = 


1.2 


M 


10 

I 


Z ^ 

D ^ ir“^ 


B 


(2.19) 


such that < I, 

i — e — t 


where 1,!^ = cracking moment calculated a.s 

M = 

yt 

M = moment on the section considered. 
Z = the lever arm 


( 2 . 20 ) 


f^“ modulus of rupture of concrete which is related 
to the characteristic cube strength f^^ aS 


fr = O-'^fou 


( 2 . 21 ) 


y.j. = distance between the centroidal aixis, of the gross 
section, to the extreme tension fibre and can be 
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computed as 

= (I> - (2.22) 

where is computed from equations (2,15) or (2.17) depen- 
ding on whether the depth of neutral axis is greater or less 
than the depth of flange, llien the short- terra, deflection is 
obtained using equations (2,9) and (2.10) replacing the 
moment of inertia (l) "by the effective moment of inertia 
(Ig) and the length (1) by the span (1^) of the beam. 

2. 3. 2, 1(b) lONG-TBHi DSFIECTIOh 

The deflection due to shrinkage and creep of concrete 
are considered in the long-teim deflection. The detailed 
analysis of these deflections are given in the article( 3. 3. 2. 1b) , 

2, 3.2.2 CRICKS 

Cracks in concrete in the tension zone of an ordinary 
reinforced structure arc unavoidable. But in case of water 
retaining structures the design is done in such a way that there 
will not be any visible cracks. The v/idth . of cracks can 
not be obtained accurately, but some empirical formulae are 
given to caLculato the v/idth of cracks (10) , detailed 

procedure for the calculation of the crack widths and their 
limitations are given in the next chapter. 
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CIUPTER III 

MiliYSIS OZf REINFORCED CONCREIE TEE BEM- LIMIT STATE APPROACH 

3.1 GBITERAL 

A structure car be built as a framed structure or 
a walled stiucture, A framed structure has beams, columns 
and slabs as the primary design elements whereas a walled 
structure has walls, beams and slabs. The most common form 
of a beam is the Tee beam. A beam constructed monolithic 
with the slab o\’-er it will function as a Tee beam. It has 
increased area in compressive zone due to the monolithic 
action of the slab, ill the codes of practice permit the 
design by the traditional and well established elastic 
theory or the ultimate load theory. Of late, the designs 
are based on the limit state theory. A design of the most 
commonly used structural element is herein presented, A 
brief summary of the properties of the materials is also 
given, 

3.2 MiTEPJAL PROPERTIES 

Concrete and steel are used as the materials of 
construction. The most widely adopted mix of concrete for 
the Tee beam and slab floors is of grade 1^0 which has a 
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characteristic cube strength (f^^) of 20N/inni^. Hot -rolled 
mild steel with an ultimate yield stress of 250 N/mm^ is 
considered in the design. Typical idealised stress- strain 
curves for concrete and steel are shown in Fige. (3.1) 
and (3.2). 

3.2.1 CONCRETE 

A rectangular parabolic stress block Fig. (3.3c) 
is used at the limit state in flexure for concrete in 
coii5)ression. The strain ’s^' in concrete at the junction 
of the parabola and rectangle is given by 

= ( 3 . 1 ) 

Tile ultimate strain s„„ in conorete is dependent on tlie 

ou. 

grade of concrete used and is given by 

®cu = (4-fcA5) X 10“5 (3.2) 

where f_ is the cylinder ( 300 mm dia and 600 mm long ) 
crushing strength. The typical relation between the cylin- 
der crushing strength and the characteristic cube strength 
is 

to = °-'78 ^eu 

The Young’s modulus (E^) of concrete is assumed to be 
independent of the stress level for a particular grade but 
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is dependent on the density (^ ) of concrete. It is 
calculated as 

Eo = 9.9<^2 (3.4) 

The long-term modulus is about half of the value given by 
the equation (3.4). 

3.2.2 STSSL 

A trilinear stress- strain relation is assumed in 
case of hot- rolled mild steel as shown in Fig. ( 3.2 ). Two 
different yield levels are assumed with the first yield 
occuring at a stress level of A and the second at B as 
indicated. Two different stress-strain relations are 
used wherein the second yield point B' in compression is 
at a lov7er level than that of tension value. The Young* s 
modulus of steel (E_) is assumed constant and same in 

o 

5 2 

compression and tension vdiich is equal to 2,0x10^ N/mm , 
3.3 STRUG TURAD ARA&ISIS 

3.3.1 LIMIT STATE OF COLLAPSE 

3. 3.1.1 BASIC ASSUMPTIONS 

The ultimate load capacity of the beam is computed 
based upon the following assumptions. 
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(1,) Plane sections remain plane even after "bending* 

(2) ITaere is a perfect "bond "between concrete and steel. 

(3) T"lie stress distribution in the concrete in flexural 
compression is given by the stress-strain curve 
for concrete in uniaxial compression. 

(4) ilhe concrete is said to fail when the farthest 
fibre is strained to an ultimate value recorded in 
a direct conipression test, 

(5) The stress in the reinforcement is given by the 
appropriate stress-strain curve obtained in direct 
tension test. 

(6) The tensile strength of concrete is ignored except 
that the stiffening effect of concrete in tension 
is considered while calculating the deflections of 
the beam, 

3.3. 1.2 PLSIURE 

Consider a typical Tee beam section shown in 
Fig. ( 3, 3 )• For design purposes it is convenient to 
assume a stress block shape which is as simple as possible. 
Many different sh^es have been suggested from time to 
time. Some of them which are in use at present are shown, 
in Fig* ( 3. 3 )• most widely used stress bloc"k: is the 

rectangular stress block proposed by C.S, Ihitney in 1937. 
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In the present thesis the rectangular parabolic 
stress block shown in i'ig, (3'. 52) is adopted. The stress 
block factors (K-j, K 2 and K^) are given by empirical 
formulae given in eq.uations ( 2.1 ) and ( 2.2 ). The total 
compressive force is 


where is the total compressive force assuming an 
ordinary rectangular cross-section of breadth equal to the 
width of the flange (B) and P^g is the compressive force 
under the flange projections, shown in Pig. (3.4a) 


^c1 = S (5.6) 

^c2 ” ®c Sc2 (5.7) 

Pq 2 < ^-1^3 fc (B-'b) (i-df) (3.8) 


where s ^2 i® strain in concrete at the soffit of the 
flange at failure, d is the effective depth of the section 
and d^ is the depth of flange. 


3.3.1.2(a) Analysis of a beam is a trial and error 
process. Recommended step by step procedure is as 
follows: 

(1) Determine the total compressive force P^. 

(2) Determine the total tensile force 

^s “ ^yt -^t ■" Asc ^yc 


(3.9) 
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vviiere Ag^ and are the tensile and con^jressive 
reinforcement respectively. 

Determine the depth of neutral axis (d^^^) 

^n1 = ^s/^o (3-10) 

1 

If (3^^ < t "felie beam is under- reinforced. 

Then go to step 5." Go to step 6 if the beam is 

/ "I 

over-reinforced (d^^ > d) , The section is 
balanced if ~ d. 

If the actual stress (f_) in steel is greater than 

o 

the first yield stress (fy^) j estimate the nev/ 
depth of neutral axis 


'^11 

Then go to step 7. 

If fg is less than ^yt- calculate f^ from the 
strain in steel (s_) 

o 

Ss/®cu = 


n 


(3.12) 


where s^„ is the ultimate strain in concrete and 

CIa 


®s ^ *^n 

where K= f^, b/A 


at ®s 


(3.13) 
(3.H) 

On substituting equations (3.13) and (5.14) in 
equation(5. 12) , an equation in terns of d^ is 
obtained from which the new approximation for 
d^ is obtained. 


( 3 . 15 ) 
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(7) Oheclc if forces and "balance; if not 

InJ = ^V^'o) <^n2 (5.16) 

(i) when under-reinforced, see if the forces 
balance for d^^. 

(ii) If over-reinforced, d^^ over corrects, hence 

try the mean of dj^g ^n3* 

(8) If forces still do not balance, estimate subse- 
quent d^ by interpolation technique. 

(g) Once the forces balance, the depth of neutral 

axis is known from which the mcmient of resistance 
of the section is determined. 


3.3.1.2(b) MEhT OP RESiS'I'ilTCB 


Refering to Pig. ( 3 . 4 b) "the total compressive force 
is divided into two with P^^ being the resultant force of 
the rectangular stress block and P ^.2 being the resultant 
force of the parabolic portion of the stress block. The 
factor (k.,) for the depth of the parabolic stress block 

is calculated. 


where f^^^^ is in h/ram 

a - V C5.18) 

Hence the depth of the rectangular stress block (dg) is 



3-11 


^2 ' " *^1 (3.19) 

If the depth of neutral axis is greater than 
the thickness of flange stress hlock given by simplified 
unified code ( 1 ) is used for the sake of simplicity and 
is given in article 3. 3.1; 2(c). 

Hence for a unit mdth of beam 

^c1 = 0*^5 f^^ d2 (3.20) 

^'c2 = -f- ^ 0-45' f^^ d^ (3.21) 

The lever arms Z^, Zg to both the forces are 
detemined from the geometry 

Z^=(d-d2/2) (3.22) 

Z2 = (a - d2 - d^) (3.23) 

Total compressive force (f^) is determined as 

]?e=r,i + re2 (5.24) 

Kno^ving the values of and Z2 the moment of resistance 
is computed 

= ^01 h + ^02 ^2 ( 5 - 25 ) 

3.3.1.2(c) Is the depth of neutral axis is more than the 

thickness of the slab, the effective area of concrete in 
compression is the full flange and the portion of the web 
above the neutral axis. The rectangular parabolic stress 
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"block poses a pro’blem of geomefry in tlie calculation of 
tile total compressive force and the le’^erajnn. So, a 
simpler and fairly accurate rcctnagular stress block ( 1 ) 
is used . Fl0' (3.3cO 

Tlie total compressive force P is divided into 
two forces, being the force in the flange and ^^2 
being the force in the portion of the web. Then 

^0 = *01 + *'o2 (3-26) 

where 

I'd = 0.4 B d^ (5.27) 

^c2 = ^cu ^ ^V^f^ 

With these values of compressive forces the depth of 
neutral axis is obtained in the same way as described 
in article 3.5.1.2(a) 

" '^n “ V2 (5.29) 

^2 = '*11 - ( ^ 

Moment of resistance is calculated by using 
equation (3.25). 

5. 3.1.3 SHEjlH 

Unlike the elastic approach wherein it is assumed 
that uniform shear stresses act on vertical sections of 
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"blie "beam* I'or tilie G''acked sec1;loii "fhe si'tu.a'tion is 
mucli more complicated. The shear is resisted hy 

(1) Uniform shear stresses on the compressive zone, 

(2) agfregate interlock along the cracks. 

(3) dowel action in the bars where the concrete 
between the cracks transmits shear force to 
the bars. 


In view of these problems the shear stress calculated by- 
simple formula similar to that given by elastic theory is 
retained 


V = 7/bd 


(3.31) 


The area of steel to resist these shear stresses 
is obtained 

_b VC. 

_ 


0.87 f 


(3.32) 


y 


The spacings (s^) of links is such that every potential 
crack is crossed by at least one link. To ensure this, 
the spacing is limited to 0.75 d in the direction of the 
span. 


Shear stresses due to torsion are out of scope of 
this thesis as the manber is not subjected to twist. 

3. 3. 1.4 BOlU) 


Bond is the frictional grip between steel and 
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concrete. Without proper bond, there is no composite material 
like reinforced concrete and the no slip condition can not be 
satisfied. Main emphasis is given to the local bond. The 
local bond is calculated in the same way as is in case of 
elastic approach. 



V 

Z So 


(5,53) 


where f-j^^ is the actual local bond stress, 

Z is the lever arm and 

So is the sum of the perimeters of all the reinforcing 
bars. 


5.3,2 LIMIT STAT3S 0? SERVICSABILITY 


3. 3.2.1 LEPLlilCTIOil 


'The limit state of deflection is concerned with 
probable maximum loading under normal working conditions 
during the life time of the structure. Hence the charac- 
teristic loads .are directly used with partial safety factors 
being 1.0 for dead load as well as superimposed load. The 
partial safety factor for materials is also 1.0, 

The dcfl3Ction of a structure has two components, 
the short-term deflection and the long-term deflection as 
defined earlier in Chapter II, The analysis for deflections 
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is based on the elastic methods and the conTenient v/ay is 
to find the curvature (<^) and then obtain the deflection(y) . 

^ (3.54) 

SI 

where M is the moment due to external loads and 

El is the flexural rigidity of the beam. 

y = k 0 (5.35) 

where k is a constant depending on the boundary conditions 
of the beam and l-j^ is the span of the beam. 

She foregoing equation (3.35) for elastic deflection 
takes no account of variation in section properties along 
the span. The support and mid- span sections, for example, 
usually contain varying amount of reinforcaaent. However, 
the use of values for these terms based on the section at 
mid-span has been found to be satisfactory ( 10 ) for both 
simple and continuous beams. 

3.3.2.1(a) SHOlS-IEKvI DEFLECTIOHS 

Elastic theory asstimes linear behaviour between 
stress and strain even at hi^ loads. The concrete in 
reinforced members hovraver, cracks at relatively low tensile 
stresses. As a result the load-deflection cu3rye tends to 
be bilinear. The short-term deflection therefore should 
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be determined according to whether or not the member is 
cracked or uncracked. But with little sacrifice in 
accuracy the deflections are measured assuming a cracked 
section. The effect of stiffening of concrete in tension 
zone is considered as given earlier in the preceding 
chapter. 

5. 3. 2, 1(b) LONG-TERM DEPISCTIONS 

The cause of long-tem deflections, as described 
earlier, are mainly due to the effects of shrinkage and 
creep in the structure. Actual determination of the long- 
term deflections involves separate considerations of the 
actual shrinkage and creep coefficients of the concrete. 
However ,jf-thee:q)ec ted shrinkage and creep strains are not 
excessive, then the total long-term deflection is assumed 
to be the sum of the short-term deflection due to all loads 
and the additional deflection due to creep and shrinkage. 

The deflections (y^) due to creep and shrinkage(yg) 
are calculate as 

Jg = P Tgg 1^ (5.36) 

v/here p is a constant depending on the support conditions, 
^cs 

v- _ ®CS 

‘cs 


(3.37) 
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where is creep ultimate strain taken as 0.0003 

-C P - p' 

S =0,72 < I1.0 

fp 

for 0,25 _< p < 1.0 (3.38) 

= 0.65 < 0 

Tp " 

for P > 1.0 (3.39) 

where p and p^ are the percentages of tensile steel and 
compressive steel respectively 


p = 

A3t 

hd 

(3.40) 

p' = 

,^00 Isc 

(3.41) 

hd 

^c = 

1 

0 

(3.42) 


where y^^ is the creep deflection associated v,ath the 
initial deflection due to permanent loads. The Young's 
modulus of elasticity is 


See = 

where 0 is the creep coefficient of concrete 

y^ is creep deflection due to permanent loads. 

3.3.3 IIIIIT STATE OE lOGAE, DAI.IAGE 


The effect of local damage due to cracks in the 
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■beams is elaborated in this v/ork. The cracks in the 
reinforced concrete beams are random and the variations 
in their width is very large. There is not much of research work 
carried out to determine, accurately the widths of cracks 
in a beam. The width of cracks are determined by en^sirical 
formulae (10 ) • ®ie surface crack width in general is 


3 S[Q 


^cr " 




(3.44) 


where = "the distance from the nearest reinforcing bar 

to the point where the crack width is going 
to be calculated as 


s = strain which is given by 
m 


= s^ - 
m ' 


1,2 b^I)(a’-d^) X 


10 


,-3 


(3.45) 




a' = the distance from the compressive face to the 
point of crack 

s. = strain at level considered ignoring concrete 
intension zone 

b. = width of the section at the level of centre of 
X 

gravity of tension steel 

e . = minimum of all values of 
min 

These oraok mdths are limited to marimuB values as 
given In the following chapter. 
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CHAPTER IV 


POaiI]I..iTIOH OP OPTIMIZATION PROBIBM 


4.1 INTRODUCTION 

Th.e invert '^ion of the digital computers has 
enabled engineers solve complex problems by numerical 
methods. Of late, various operations in structural 
design have been subject to extensive research, the most 
significant being the optimization of structures, A 
design does not mean simple proportioning of the relevent 
components of the structure but should give a structure 
which is optimum and efficient. Optimization of the 
most common element in structures namely the Tee beam 
has been attempted herein. There are ' number of 
mathematical methods which enable an engineer to optimize, 

4,2 G-MERAD PORIIULAI'ION AMD METHODS OP SOLUTIONS 

A structural system can be described by a set 
of quantities which are viewed as variables during the 
design process. Tbe quantities that are fixed at the 
outset are called preassig^d parameters, while those 
quanities that are arrived^during the process of design . 
are called design variables. 
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The oojective function in a structural design is 
the basis for choice amongst all the alternative designs, as 
the design has infinite number of solutions . The opti- 
mization of a total structure as a whole bg s bearing on 
the optimization of each element of the structure* It 
is difficult to select the objective function vdiich effecti- 
vely relates the element to the whole structure* So, one 
of the guiding factors to select the objective function 
may be the one which is an important design property that 
can possibly be quantified,, 

A general optimization problem can be stated as 
find D , 


that minimizes an objective function 



F(D) 

(4.1) 

such 

that l;| 5 -(jD) = 0 for fc= 1,2,,**,l5: 

(4.2) 

and 

g.(D) < 0 for j = k+1, K+2,,..,m 

J 

(4.3) 


The design vector I) is a point in n-dimensional design 

space, function l^^-CB) denotes the equality constraints 

and g-(D) the inequality constraints* The objective 
J 

function is I (D)* 

Ifeny classical methods and numerical methods 
ane available to seelc solution of the problem stated in 
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equations (4.1), (4.2) and (4.3). The difficulty in 
solving the problem by classical methods necessitates 
the application of numerical techniques. There are 
algorithms like linear programming, non-linear programming, 
geometric programming, dynamic programming etc. One of 
the methods to solve such a problem is a linear progra- 
mming method wherein the equations (4.1), (4.2) and (4.3) 
are suitably replaced by linear approximation by Taylor 
series expansion about the point D, neglecting all the 
quadratic and higher order terms. This alternative 
enables the application of well established linear progra- 
mming algorithms to solve the given basic problem. 

There are other methods #iich can handle the 
problem as given in the preceding paragraph. The methods 
like penalty function formulations are widely used methods 
for nonlinear problems where the objective function and 
the constraints are nonlinear functions of the design 
variables, 

4.2.1 PMAliTY FUNCTION METHODS 

Of the methods of optimization of a nonlinear 
constrained problem, penalty function methods are the 
most widely used. Penalty function method transforms the 
basic constrained optimisation probla. into alternative 
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formulations such that numerical solutions are sought hy 
solving a sequence of unconstrained minimization problems. 
The general form of formulation of an unconstrained minimi- 
zatioii problem can be given as 

, -* m r” -► _► 

0k = f CD) + S lk(D) , g.(D) (4.4) 

i=1 L 3 J 

where Grj_ is some function of equality and inequality 
constraints I^Cd) and gj(D) respectively and r^ is the 
penalty panameter. If the unconstrained minimization of 
is repeated for a sequence of values of , 2 , 3 ...) , 

the solution may be brought to converge to that of the 
original problem stated by equations (4.1) , (4.2) and 
(4.3). This is the reason why penalty function methods 
are also known as sequential unconstrained minimization 
techniques (SUIH). 

The algorithms available to adopt the technique 
of SUMT are many in number. They can be classified into 
board catagorics as direct search methods and descent 
methods. The direct search methods do not require the 
partial derivatives of the given function to be optimized, 
while the descent methods need the gradients, for the 
same reason the direct search methods are also called 
nongradient methods. ComparatlTely dasoent methods are 
more efficient if the function consists of a large number 



4-5 


of variables* The dcscont method called ’Variable Metric 
method' or the ’Davidon fl etcher Powell’ (DPP) method 
( 14 ) is used in the present optimization. The DPP 

method is used along with cubic inteipolation method for 


one- dimens tonal minimization, to obtain the step length 
^opt determine the new design vector. 



Di-. 


\pt ®i 


(4.5) 


The interior penalty method is used 
Jn the penalty methods to solve the constrained 
problem. The main advantage of this being that the vector 
obtained after every SUIT is in the feasible domain of 
design ^ace. But it is necessary that this method has to 
be supplied with an initial vector vMch is in feasible 
domain of design space. In case of engineering problems 
this does not pose a great problem as it is not much very 
difficult to supply a Dq is in the feasible domain. 

There is also a method to obtain an initial feasible point 
using the same penalty function but it is out of scope of 
the present work. The fornulation in case of the present 


problem is 





1 

i^D) 


( 4 . 6 ) 


as there are no equality constraints. 



4.3 IOK'HILaTIOI'T OP 2HE PROBLM 


In many practical problems, the design variables 
Can not be choosen arbitrarily; rather they have to satisfy 
certain specified functional and other requirements.' The 
restrictions that must be satisfied in order to produce an 
acceptable hesign are collectively called design constraints. 
The constraints which represent limitations on the behaviour 
or performance of the system are called behaviour or func- 
tional constraints. The constraints which represent physical 
limitations on the variables are called geometric or side 
constraints, '-Tith the foregoing in view the constraints 
are divided into two groups and are given as follows. 

4.3.1 BSHA'/IOUR COISTRillTS 

The constraints that are to be considered in 
this problem are listed in the follovdng. 

( a) The limit state of moment (IJ^) in the slab due 

to the physical influence should be less than or 
equal to the moment of resistance of the 

slab section so that the slab is safe against 
flexural failure. 
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("b) The T'ee boom should have sufficient floxural 

strength to resist the limit state moment (Mg) 
caused due to the external loads 

% i "Vb (4-8) 


(c) 


(<i) 


(e) 


The total shear force should safely he resisted 
by the shear capacity of the beam 


< 

B - 


all 


(4.9) 


The beam should have sufficient stiffness 
associated with the flexural strength so as 
to be functionally sound and should limit the 
deflection (y^,y 2 ) to a maximum of span (li^)/250 
in case of total load and to span (l-j^)/350 in 
Case there is the dead load plus the permanent 
load. The deflection is the total of all due 
to shrinkage, creep, teirperature including the 
initial elastic deflection. 

yi i ^^all'l 

where (yaii)-| = 

^2 - ^-^ 311^2 

where ( 7311)2 = 20 mm 

whichever is less 
It is also necessary that the tensile craclcs 
in concrete in the tension zone are limited to 


( 4 . 10 ) 

(4.11) 

(4.12) 

(4.13) 
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(f) 


(s) 


accoant for the functional utility of the structure 
sUKi to safeguard against corrosion of steel and 
spoiling Of concrete, Ihe m^-^xinium surface crach 
width (^) is limited to 0.3 mm and the sui-face 
crack width (wj just belov. the reinforcing hor 
to 0.4- dg 

( 4 . 1 ^-) 

( 4 . 15 ) 


W 

m 


< 0.3 


W3 < 0.4 dg 


fo avoid sudden failure of concrete and to inpn.rt 

ductility to the structure depth of nouural 

\ - 1 ■ X +-hp most to the balanced 

a:ric (d^^) is liraited at the 11 ^ 

or a limiting depth of meutral axis (dj^p) 

- 0.0715 (4.16) 


^nl “ 


250 


^cu 


d < d T 

n 1 nl 


( 4 . 17 ) 


, ,1 be less than or 

The local bond stress shoulu 

T 4 . -rn V-. -hnnd stress to have no 

equal to the allowable bonu 

slip condition 

( 4 . 18 ) 


^b < 


^ba 


4 . 3.2 

(h) 


GBOIIETRIC COhSlEilHTS 


.-nns the code spefifics 
For all practical condition 

slab, which in turn is 


a limit on the depth of 


the depth of the flange, 
^f 


to 73 mm. 


75 < d. 


( 4 . 19 ) 



To place the reinforcement effectively the width 
of the rih of the Tee beam has to he limited to 
150 mm Or so, 

150 < h (4.20) 

Stiffness limitations on the slab as specified 
hy code . (11 ) limits the span (Ig) to depth (d^) 
ratio of the slab to 26. 

Ig/df < 26 (4.21) 

The span of the beam (l-j^) is arbitrarily ' limited 
to 10,000 mm as it is 'uneconomical to have larger 
spans, 

1-b i (4.22) 

The code specifies the limitation on span (1^) to 
depth (d) ratio to control the deflection of the 

beam, rigorous analysis is done the span to 

depth ra^tio is arbitrarily limited to a lesser 
value than that of the code to make this constraint 

an inactive constraint, 

1 ^/d < Df (4.23) 

where for cantilever beam 

= 15 for simply supported beam 
= 20 for continuous beam (4.24) 

To take care of temperature effects and to control 
the cracks there is a minimum amount of longitudinal 
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(Ag-jjv^) and reinforcement (A^^) to be provided. 


^stm 

^t 

(4.25) 

- 

^sv 

(4.26) 

^stm 

=0.85 bd/fy 

(4.27) 

A 

^svn 

= 0.002 b c^ 

(4.28) 


is the spacing oi shear links, 

(n) The spacing of links should be so that every 
potential crack is crossed by a link. To ensure 
this a maximum spacing of 0,75d is stipulated 

c^_< 0.75 d (4.29) 

(o) In no case a design variables are less than zero 
as they are all physical quantities 

< 0 (4.30) 

i=1,2,3, ... ,11 

4.3.3 THE DESIC-N VECTOR iHD THE COST EUIfCTIOiT 
The design vector is 

D = (xji=1,2,3, ,n) (4.31) 

where n is, the number of design variables, equal to 

11 in the present problem. 

Refering to Eig.( 3.3) the components of the 
design vector are 
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= B 

Xg = "b 

X3 = 

^10" °st 

r 


x^= d 


= A. 


^ ^sc 

^ " Asv 

^6"" Xg = Ag^(slab) 


^5= D 


(4.32) 

where all terms are as defined earlier. 

The foregoing vector is to be determined subject to the 
constraints given in articles 4.3.1 and 4.3.2, optimising 
the cost function or the so called objective function. 

Total cost = 3?(D) = j^Bd^ + b (Ii-d^)J l-j^ 

^r^ ^b ^1 

+ (li^ 2 d — 2 dj) R2 Bg ^c^^b (4.33) 
Uondimensionalizing it 

Cost/Gg = B d^+ b(D-dj) l-,^ 

^-^sc+-^sf^^v^ ^b ^ 

+ (l.t,+2]>2d^) Rg Ig/l^ (4.34) 

where ^s'^'c (4.35) 

IT is the density of steel and 0 is the 

s 

unit cost of steel ils/Newton. 


R2 


G lx, 
w b 

”g^ 


G is the unit cost of v^ood Rs/mm, 
w 


(4.56) 
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4.3.4 3'0RI.Iin.ATI01T-SMDAJ?D FOEM 


Find vector F i=1 , 2, 3, . . . . , 1 1) to 
F(D) = |b + d (]>-d^)] 

+ (i.gQ + Ag^ + Agy) 1-^ R-| 

+ (1^ + 2F - 2d£) Eg lg/1^ 


subject to 


gl Mg / i^g - 1 < 

gg = Mb / - 1 < 

g3 = vb / y^2r ^ - 

S 4 " y-i / ~ 

% = ^2 /(yall^2“^ - 
gg = / 0.3 -1 i 

g^ = Wg / 0.4 dg-1 < 

§8 <in / - 

% = "^b / ^ba ' ^ 

gio= ^ 

150 / b - 1 < 

S 12 = ^s / - 

g^3 = lb / 10000-^ < 

gl4 = lb / V " ^ “ 

S'i5 " "^stn/ - 

S 16 “ -^m/ '^r “ ^ ~ 
/0.75 d-1 < 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


minimize 


(4.37) 


0 
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for i=1,2,3,*...,11 

and 3=18,19,.. .,27 (4.38) 

- The problem stated in the foregoing equations 
(4.37) and (4.38) is solved by interior penalty function 
method. The results and discussion are presented in the 
following chapters. 
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cost of steel is assuj-ned to te Rs, 3.50 , 4.00 or 4,50 per 
kilograia. This accounts for the fluctuations in the costs 
of these vital materials of constiruction. The variation 
in the cost of form work is not taken into account as the 
hype of material chooseii for form work differs greatly 
from place to place and usually the cheapest material is 
used, A constant value of Rs,25 per square meter is assumed 
as the cost of form work. It is fairly easy to incorporate 
the variation in the cost of form work if sufficient data 
is available, A set of nine problems are solved for given spans 
of the beam and the slab. The cost ratio (OR) versus the 
total cost, per millimeter span of the beam, for different 
span ratios are plottod as in Figs, 5.1 through 5.7. The 
cost ratio is calculated as : 



where j is the density of steel, C and are the costs 
of steel and concrete respectively. It can be observed 
that different curves result in for different costs of steel. 

The sp^ ratio between the beam and the slab is 

given by 

Span ratio = li/lg (5,2) 

Table 5.1 is drawn giving the areas of steel for a 
of slab (3000 mm), with different spans of the 


given span 
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beam. The costs of concrete and steel of Rs. 4 OO per cubic 
meter and Rs» 4 per kilogram respectivel 7 , are adopted for 
the optimization of the area of steel in the Table 5.1. 
Table 5.2 gives the areas of steel for a single span ratio 
and different cost ratios. 


TABI.E 

5.1 

: AREiS OF STESI FOR DIFFERMT SPANS OF THE 

BEAM FOR 3000 MM SPAN OF THE SLAB 

S.No. 

Span 

the 

(mm) 

of 

beam 

Area of 
tension 
steel 
(mm 2) 

Area of 
compressive 
steel 
(mm^) 

Area of shear 
reinforcement 

(mm^) 

1 

3600 


7995.57 

0,0 

990.17 

2 

4000 


7995.57 

0.0 

990.17 

3 

4500 


7995.57 

0.0 

990.17 

4 

5000 


7995.57 

0.0 

949.00 

5 

6000 


7996.13 

0.0 

1092.68 

6 

7000 


7996.13 

0.0 

1092.68 

7 

8000 


7996.13 

0.0 

1092.68 

8 

9000 


7996.13 

0.0 

1092.68 
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TAELE 5.2 : AREAS OP STEEL FOR DIFFERHITT COST RATIOS 

THE SPAR RATIO IS 2.5 


S.No. 

Cost of 
concrete 

(Rs/m^) 

Cost of 
Steel 

(Rs/kg) 

Area of 
tension 
steel 
(mm2) 

Area of 
c empress ive 
steel 
(mm^) 

Area of 
shear 

reinforcement 

(ram^) 

1 

400.0 

5.50 

7989.75 

0.0 

774.82 

2 

500.0 

5.50 

7989.75 

0.0 

774.82 

3 

600.0 

5.50 

7989.75 

0.0 

774.82 

4 

400.0 

4.00 

7989.75 

0.0 

774.82 

5 

500.0 

4.00 

7989.75 

0.0 

774.82 

6 

600.0 

4.00 

7989.75 

0.0 

774.82 

7 

400,0 

4.50 

7989.75 

0.0 

774.82 

8 

500.0 

4.50 

7989.75 

0.0 

774.82 

9 

600,0 

4.50 

7989.75 

0.0 

774.82 


It is interesting to note that the area of steel is 
almost the same in all the cases. This can he visualised 
hy studying the results presented in Table 5.1 and Table 5.2. 

The initial feasible point given and the optimal 
design vector obtained are sho^wn in Table 5.5 for a typical 
example with span ratio of 2.855. It is observed that the 
initial feasible point is close to the optimal point. Of 
this reason, the problem has converged very fast and at the 
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mos'fc ilj "took tlaree cycles of SUMT witli the initial value 
of the penalty parameter = 1,0, 

TABLE 5.3 : T4BLE SHOWING INITIAL FEASIBLE POINT 

AND THE OPTIEiAL DESIGN VECTOR (SAPN RATI 0=2 ."83 5) 


Design 

Vector 


^3 

^4 

®5 

^6 

^10 

^11 


Initial feasible 
point - 
■^0 

Optimal design 
vector 

^opt 

0.0 mm 

1124.96 mm 

200,0 mm 

1 51 » 71 mm 

200.0 mm 

115.54 mm 

950.0 mm 

1219.45 mm 

990,0 mm 

1259. §-5 mm 

8000,0 mm^ 

7989.86 m? 

2 

0,0 mm 

2 

0.0 mn 

150,0 mm^ 

722,82 mm^ 

2 

20,0 mm 

683.12 mm^ 

50,0 ram 

668.20 nm 

200.0 ram 

461.32 mm 


The same initial point is choosen for all span ratios. But 
the obtained optimum point is successively used as the initial 

point for other cost ratios. 
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5.3 CONCLUSIONS 

(1) The hierarchy of the values of the critically satisfied 
constraints in general is observed to be in the oi^er as 
follows: 

( i) the limit of depth of natural axis, tending the 
design to be a balanced design, 

(ii) the limit on the span to depth ratio of slab, 

(iii) the limit on the depth of the slab, 

( iv) the limit on the minimum width of the rib, 

(v) the limit on the spacing of the shear reinforcement, 

(vi) the limit on the span to depth ratio of the beam, 

(vii) the moment capacity of the slab., 

(viii) the deflection of beam., 

(ix) the limit on the minimum area of shear reinforcaaent, 

(x) the moment capacity of the beam, 

(xi) the limit on the deflection y.| of the beam, 

(xii) the maximum crack width, 

(xiii) the average crack width and 

(xiv) the shear capacity of the beam. 

The constraints specifying the nonnegativity of the variables 
are neglected. It is also observed that almost all the 
examples followed the same hierarchy. Jill tho constraints 
that are critically satisfied are the geometric constraints. 
This infers that the beam is sound in its behaviour. Hence, 
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"tlie codes can "be relexed in constraining the geometry of 
the structural elements in ordinary buildings vAiere there 
is not heavy loading. 

(2) As the computed deflections are less than the 
permissible deflections the span to depth ratio in the 
codes Can be very well relaxed as the limitation on the 
span to depth ratio on the beam is arbitrarily fixed in the 
examples under consideration. It can also be concluded that 
the optimal beams have high flexural rigidity. 

(5) As there is no limit prescribed on the maximum depth 
of the beam the optimal solutions will have ribs of large 
depths. Now-a-days architects tend to adopt diaphram tjrpe 
ribs in case of beams. It is so observed that all the solutions 
giving optimal design for cost invariably have given diaphram 
ribs. But a rigorous analysis is to be done for the lateral 
buckling and the effect of shear on the deep beam. 

(4) Yariation in costs is assumed for concrete and steel 

and the curves indicated that the total cost as a function 
of the span of the beam is less when the cost ratio is 
larger. Bor given costs of steel and concrete the total cost 
of the beam can easily be cemputed from the graphs. As an 
examiple let us assuime the market price of steel to be Rs,4.0 
per kilogram and that of concrete to be Rs. 4^X).0 per cubic 
meter. Assuming the density of steel to be 7.85 the 
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cost ratio OR is 

0,0078 X 40 

CR= — *» 7 eo-,o (5.5) 

400 X 10”^ 

where the quantities are expressed in terms of Newtons and 
millimeters in the equation (5.5). 

Then, the total cost of structure Can he obtained for a given 
span ratio for the optimal problem, Assume a span ratio 
of 2,5, Hence the total cost of the structure frOT Rig.5.5 
Can he obtained as 0.22 per millimeter span of the beam. 
Knowing the span of the beam in millimeters say 5000 mm, 
we can calculate the total cost as 

Total cost = 0.22 x 5000 Rs. 

= Rs, 1100.0 (5.6) 

(5) Obserrtng the hierarchy of the values of constraints, 

the more critical heing the geometric constraints, it can he 
oonoluded that the codes are still conservative and a modifi- 
cation can alWS he incorporated. Bit it should he home 
in mind that a complete check for slenderness is yet to he 


done 



Total Cost(Rs/mm span of beam) 
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CHAPTER yi 


OViE Vim 


6.1 GMERAL 

The traditional methods that are in use are not 
rational regarding the actual factor of safety. It is 
a well known fact that the actual feet or of safety in 
case of the well established elastic method of design 
varies with the loading and structure. The plastic methods 
developed later have considered this factor and have 
suggested load factors to compute the design loads. Even 
then these plastic methods have not reflected any light 
about the structural behaviour at working loads. But it 
is absolutely necessary for a structure to be rationally 
designed and checked for all serviceability conditions. 

So a method called the "limit state method" is the 
culmination of all the methods evolved hitherto. It 
incorporates the load factor design with due account of 
analysis of the structure for working loads. This method 
of design tes been taken up and a Tee beam element is 
optimized. This is a step in the right direction of 
optimizing large systems. The present method of analysis 





can further be improved. 

6.2 SUGGESTIOHS POE PUTUHE ¥(HK ; 

1 - The same type of work nay be repeated varying 

the grade of concrete and the nature of steel 
reinforcement to geu exhaustive results for 
various types of reinforcement and grades. 

2. The rectangular parabolic stress block can further 
be improved by incorporating Hogenstad's (1) 
stress block. 

3. The member can be analysed as a continuous 
member, with due consideration of the different 
moments of inertia at the midspan and supports. 

4. The redistribution of the bending moments and 
different flexural rigidities have to be incorporated 
in case of continuous members. 

5- The bilinear nature of the moment curvature 

relation should be accounted for computing the 
deflections of the structure with uncracked and 
cracked section. 

6. The total structirral system may be optimized 

along v/ith the elements of it and study of the 



cost in both the cases should be studied. 

A detailed analysis with different ways of 
reinforcement in the estimation of crack widths, 
can be carried out. 

The study may be extended with the application 
of probability theories . 

Optimization may be done with different methods 
of programming like Geometric programming, 
Dynamic programming etc. and the results be 
compared for further insight. 



R-1 


HEPFRSTCES 

1. Barry P. Hughes; "Limit state theory for reinforced 
concrete", Pilmon Publishing, 1971. 

2. Bengt B. Broms; "Crack v/idth and spacing in R-C. 
members"; Jl. Amer. Cone. Inst., V 62, 1965, pp .1237-1257 . 

3. ACI Committee Report; "Control of cracking in concrete 
structures " 5 Jl . Amer. Cone. Inst., V 69/12, 1972, 

pp 717-754. 

4. Peter Gergley; "Distribution of reinforcement for 
crack control"; Jl. Amer. Cone. Inst. 

5. Edv?ard G. llawy; "Crack control through reinforcement 
distribution in two-way acting slabs and plates"; 

Jl. Amer. Cone. Inst.; V 69/4, Apr. 1972, pp. 217-220. 

6. Glen S. Orenstein and Edward G. llawy "Crack Viridth 
control in R.C. two-way slabs subjected to a 
uniformly distributed load"; Jl . Amer. Cone. Inst. 

V 67, 1970, pp. 57-61 . 

7. Robert H. ELvery and Mohamad Shafi; "Analysis of 
shrinkage effects bn R.C. structures"; Jl. Amer. 

Inst., V 67, 1970, pp- 45-53. 

. Dan E. Brauson; "Deflections of R.C. flejoaral 

members", Jl. Amer. Cone. Inst.; V 63/1; Jan 1966, 
pp. 637-674. 


8 



R-2 


9. Donald G. Alcock and Andinan Pauw; "Controlled 
deflection design method for E.C. beams and slabs"; 
Jl. Amer. Cone. Inst.; V 59 j May 1962, pp. 645-659 . 

10. T.J. MacGinly; "Reinforced Concrete-Design theory 
and examples"; E & E.il. Spon Ltd., 1978. 

11. CP 110 : 1975; "Code of Practice for structural use 
of concrete"; 1975. 

12. Cement and Concrete Association; "Handbook on the 
Unified Code for structural concrete" (CP 110:1972); 
1974. 

15 . IS : 456 (Draft Code) : Code of practice (Third 

revision) for reinforced concrete structures, 1978. 

14. Rao, S.S; "Optimization - theory and applications"; 
Wiley Eastern Limited, 1978. 

15. G.I.H. Rozvauy; "Optimal design of flexural systems"; 
The P egamon Pres s ; 1976. 



